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ABSTRACT: 

A  technique  for  calculating  ballistic  winds  from  satellite  infrared  spectro- 
photometer (SIRS)  radiance  data  has  been  tested.    The  experiments  with 
SIRS  data  from  both  NIMBUS  III  and  IV  show  the  technique  has  consider- 
able skill  in  computing  ballistic  winds. 
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Detachment  FAMOS 


ABSTRACT 
A  proposed  ballistic  wind  determination  technique  using  the  satellite 
infrared  spectrophotometer  (SIRS)  data  from  NIMBUS  III  and  IV  was  tested. 
The  general  approach  was  to  derive  regression  equations  from  a  dependent 
data  sample  of  nearly  coincident  subsatellite  points  and  rawinsonde  loca- 
tions .    The  validity  of  the  equations  was  checked  with  a  smaller  independ- 
ent sample  of  SIRS  radiances  from  a  later  period.    Tests  with  the  ballistic 
wind  component  normal  to  the  path  of  the  NIMBUS  III  suggested  standard 
errors  of  less  than  10  kts  can  be  achieved.    Both  zonal  and  meridional  com- 
ponents can  be  calculated  from  SIRS-B  radiance  data  which  is  grouped  in 
clusters  .     Because  the  geostrophic  ballistic  wind  is  proportional  to  a  hori- 
zontal gradient,  several  techniques  were  attempted  for  calculating  the 
gradient  at  the  central  point  of  the  irregularly  spaced  set  of  data.    With 
proper  data  handling  and  adequate  spacing  of  the  SIRS-B  soundings,  the 
accuracy  of  the  ballistic  wind  calculations  should  also  be  near  10  kts. 
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I.     INTRODUCTION 

In  an  earlier  report,  Elsberry  and  Martin  (1971)  tested  the  feasibility 
of  using  cloud-corrected  radiance  measurements  from  the  s_atellite  infra- 
red s_pectrometer  (SIRS)  instrument  onboard  NIMBUS  III  to  directly  calcu- 
late the  ballistic  density.    The  basis  for  the  statistical  technique  was 
suggested  by  K.  W.  Ruggles  CDR  (USN)  who  also  proposed  that  ballistic 
wind  components  could  be  calculated  from  the  SIRS  radiances  in  regions 
in  which  the  geostrophic  assumption  is  applicable. 

The  ballistic  wind  is  the  weighted  integral  of  the  wind  encountered 
at  the  various  levels  of  the  atmosphere  through  which  the  projectile 
travels ,  that  is: 

Vb=t-^ |    P°V(lnp)    W      (lnp)    d(lnp)  (1) 

In  p0         J      0    v  v 

where    v      is  the  ballistic  wind  through  the  depth  of  atmosphere  with  sur- 
B 

face  pressure  p„  .  and  v  and  W    are  the  actual  wind  and  ballistic  wind 
°  v 

weighting  factor  at  a  pressure  level  p.  In  general  W  is  proportional  to 
the  percent  of  the  time  a  re-entry  vehicle  remains  in  that  layer  (Finke 
1969).  Present  methods  of  evaluating  the  ballistic  wind  involve  the  use 
of  rawinsonde  and  rocketsonde  data  for  the  levels  in  the  integration  of 
Eq .  (1).  Description  of  the  proposed  ballistic  wind  determination  tech- 
nique from  satellite  measurements  is  given  in  the  following  section  (see 
also  Appendix  A  of  Elsberry  and  Martin) . 

Wright  (1971)  has  tested  three  different  schemes  for  calculating 


.  i 


ballistic  winds  from  SIRS  radiances  from  both  NIMBUS  III  and  IV.    The 
general  approach  was  to  derive  regression  equations  from  a  dependent 
data  sample  of  nearly  coincident  subsatellite  points  and  rawinsonde 
locations.    Then  the  validity  of  the  equations  was  checked  with  a  smaller 
independent  sample  of  SIRS  radiances  from  a  later  period.    Wright's  re- 
sults for  the  SIRS-A  experiment  (NIMBUS  III)  and  SIRS-B  experiment 
(NIMBUS  IV)  are  summarized  in  Sections  III  and  IV.    In  the  last  section 
some  techniques  to  improve  the  ballistic  wind  determination  from  the 
irregularly-spaced  SIRS-B  locations  are  examined. 


II.  DEVELOPMENT  OF  EQUATIONS 
Ac     FEASIBILITY  OF  THE  USE  OF  SIRS  DATA 

If  the  geostrophic  approximation  is  assumed,  Eq .   (1)  can  be  ex- 
pressed as: 

Vb  =  ~r F  P°    -9    lk  x  v„  Z(x,y,lnp)W   (lnp)d(lnp)  (2) 

v  In  p  ^  7  p  v 

0  o 

or,  rearranging  terms: 

VD  -9  Ik  xv        r  P0     Z(x,y,lnp)W     (lnp)d(lnp)  (3) 

O  0 

where  g  is  the  acceleration  due  to  gravity,  f  is  the  Coriolis  parameter, 
and  V  Z  is  the  horizontal  gradient  of  geopotential  height  on  a  constant 
pressure  surface.    Expressing  the  weighting  function  W    as  a  linear 
combination  of  known  arbitrary  functions  of  In  p,  say  F.   (In  p) ,  we 

may  write 

n 
W     (In  p)  =    L        a.  F.     (In  p)  (4) 

i  =  l      l     1 

where  a.  are  constants.    Similarly,  Z(x,  y,  In  p)  can  be  expressed  as 

n 
Z(x,  y,  In  p)  =  S     r.  H.     (s ,  y.  In  p)  (5) 

i  =1    1     1 

where  r.  are  the  constants  and  H    are  a  set  of  arbitrary  functions  of 
i  i 

In  p.    Substitution  of  Eq .   (4)  and  Eq  .   (5)  into  Eq.   (3)  produces 


n 

P 


Vd  -9  lkxV     E        a.  r.    »Po  H.  (x,y  ,lnp)F.  (lnp)d(lnp)  (6) 

NVB"flnpn  P  i    =  1  1    1  J     e     l  l 


Now  consider  the  radiative  transfer  equation  for  the  radiance  I 

i 

J.T- 

observed  by  the  i      channel  of  the  SIRS,  sensing  at  wave  number  i>  . 

i 


I.  =  I    .  r 


o1 


A,y,T(lnp),<Kv,lnp)  ]+JpB.  j~y,T(lnp)  1    d  t  {v,  lnp      d(lnp)  (7) 

0    1L  J  d(ln  p) 


Here  I      .  is  the  radiance  emitted  by  a  contaminant  such  as  the  top 
surface  of  clouds  expressed  as  a  function  of  A  (amount  of  cloud  cover), 
v,  and  T(ln  p) ,  the  temperature  structure  above  the  clouds.    The  Planck 
radiance  at  wave  number  yand  temperature  T  is  represented  by 
B.      v,  T(ln  p)    i  ,  and  r  (v  ,  In  p)  is  the  fractional  transmittance  of  the 
atmosphere  at  wave  number  i/from  pressure  level  p  to  the  SIRS.    Smith 
et  al  (1970)  have  shown  that  the  cloud  contaminant  term  can  be  evalu- 
ated by  an  iterative  procedure.    With  that  term  evaluated,  Eq.   (7)  becomes 


N.  =  P  °       B.    V,  T  (lnp)  d  T  (i/,lnp)  d(lnp)  (8) 

1       "   P°      1L  J  d    (lnp) 

where  N    is  the  radiance  corrected  for  cloud  contamination  effects ,  also 

i 

called  the  clear-column  radiance.    Since  the  only  constraints  on  F.  are 

i 

that  they  are  known  functions  of  In  p,  we  specify 

F.  =       d  t     (v,  In  p) 

1  d  In  p 

because  the  envelope  of  the  functions  F.  has  the  same  shape  as  the  ball- 
istic weight  factors.    Similarly,  since  B.  \  v  ,  T  (In  p)     is  an  increasing 


function  of  T  (lnp)  and  satisfies  the  constraints  put  on  H    in  Eq .    (5) 

i 

and  we  can  set: 


H.  =  B.    v  ,  T  (In  p) 
\        1  — 

Substituting  these  expressions  into  Eq.    (8)  yields 


N.  =  H.  F.    d(ln  p) 

l       J  p  li 


(9) 


Using  this  in  Eq  .   (6)  results  in  a  direct  expression  for  ballistic  winds 
in  terms  of  the  SIRS  radiances 


VR=. 


B 


f    lnp 


n 
Ik  x  v       L      a.  r.  N. 

Pi-I      Ul 


(10) 


where  n  is  the  number  of  SIRS  channels  used.     In  practice,  the  coeffic- 
ients a.r.  are  grouped  together  and  are  statistically  determined  through 
the  use  of  the  BIMED  02R  stepwise  multiple  regression  routine  (Dixon, 
1966).     This  produces  a  final  equation  of  the  form 


w     =-1     lkxv  .&.    b.     N.     +  C. 

lfB         j  P    Li=l      l       i  1 


(11) 


where  b.  are  the  regression  coefficients  and  C.  is  the  regression  equat- 
ion  constant.    The  constant  terms  g  and  In  p     are  absorbed  by  the  co- 
efficients ,  b.  . 


B0     BALLISTIC  D  VALUE  APPROACH 

The  geostrophic  assumption  is  explicit  in  Eq.    (11),  and  we  may  inter- 
pret the  term  within  brackets  as  a  geopotential.     Ballistic  wind  weight- 
ing factors  applied  to  the  height  deviations  (D  values)  from  standard  geo- 
potential heights  were  used  in  calculating  the  bracketed  term,  which  we 

call  the  ballistic  D  value.    This  method  of  calculating  W   then  consisted 

"  B 

of  two  steps.    First  a  regression  equation  for  the  D  value  was  calculated 
from  a  dependent  sample  of  rawinsonde-derived  D  values  and  SIRS  rad- 
iances .    All  available  data  at  mandatory  reporting  levels  between  the 
surface  and  10  mb  were  used  in  calculating  the  ballistic  D  values  from 
rawinsonde  observations,   (see  Elsberry  and  Martin  1971  for  discussion  of 
the  application  of  the  ballistic  weighting  factors  for  layers  between  man- 
datory levels  from  the  values  given  by  Finke  1969).    Above  the  termin- 
ation level  of  the  sounding  the  temperature  values  were  extrapolated  to 
1  mb  using  the  standard  January  lapse  rate  at  45  N,  and  the  hydrostatic 
equation  was  used  to  complete  the  D  value  profile.    The  second  step  was 
to  use  Eq .   (11)  with  the  derived  regression  equation  to  calculate  the  ball- 
istic wind. 


C.     GRADIENT  OF  RADIANCES  APPROACH 

Rearrangement  of  the  order  of  Eq .    (11)  yields 


V-r 


1     r   "        c.     (Ik  x  vTTN.)  +  C    1 


T         i  Hi'  2 

1*  1 


(12) 


Although  mathematically  equivalent,  the  procedure  for  evaluation  of 
Eq .   (12)  is  different  than  the  previous  approach.    In  this  case  the  re- 
gression equation  was  derived  by  correlation  of  values  of    V      with 

B 

gradients  of  the  individual  radiances  N.  as  shown  within  the  parentheses 

Thus  the  regression  coefficients  c    and  the  constant  C„  are  not  equiv- 

i  2 

alent  to  b.  and  C    in  Eq.   (11) .    A  dependent  sample  relating  the  observed 

east-west  component  (u  )  and  north-south  component  (v  )  was  obtained 

B  B 

from  Eq .  (1)  using  rawinsonde  observations  with  nearly  coincident  SIRS 
radiances.  Missing  interior  mandatory  level  winds  were  linearly  inter- 
polated. Completion  of  the  sounding  to  1  mb  was  by  extrapolating  par- 
allel to  a  mean  zonal  wind  profile  for  that  component  and  by  decreasing 
the  v  component  by  five  percent  with  each  mandatory  level  increment. 
Since  soundings  were  not  processed  unless  data  was  present  to  at  least 
100  mb,  the  reduction  toward  zero  for  the  v  component  is  probably  just- 
ifiable . 


III.  SIRS -A  EXPERIMENT 

A.   DATA  AND  DATA  HANDLING 

The  data  used  in  this  experiment  were  from  the  same  set  used  by  Els- 
berry  and  Martin  (1971).    Radiances,  corrected  for  cloud  contamination, 
were  available  for  the  eight  SIRS  channels  during  two  four-day  periods: 
December  27-30,  1969,  and  January  2-5,  1970.    The  area  of  data  cover- 
age was  the  European  continent  including  the  area  of  the  U.S.S.R.  west 
of  the  Urals,  and  the  time  period  covered  was  from  0800  GMT  to  1200 
GMT  for  each  day.    The  dependent  ballistic  D  values  and  ballistic  winds 
were  computed  from  the  1200  GMT  rawinsonde  data  from  the  same  periods 
The  rawinsonde  data  did  not  include  any  data  from  France  or  the  British 
Isles.    In  some  cases  this  was  a  limitation  due  to  the  location  of  sat- 
ellite subtracks  lying  in  this  data  sparse  zone.     In  addition,  data  was 
typically  quite  sparse  in  the  Mediterranean  region. 

The  2  7  degree  longitudinal  spacing  of  NIMBUS  III  SIRS-A  data  pre- 
cludes the  calculation  of  both  u  and  v  components  of  ballistic  winds  by 

the  geos trophic  approximation.    Instead,  V      ,  the  wind  component  nor- 

BN 

mal  to  the  satellite  subtrack  was  calculated  using  Eqs.   (11)  and  (12)  mod- 
ified as  indicated  below 


VBN  =  "f    |     [    l       b.N.  +  C.    1  ,,,„, 

f     2n    L  .*<,       i     i         1  J  (J  la) 

i=l 


VBN=  E  d,     dNi        +    C0  (12a) 


1      T    S  d.     5Ni        +    C0    1 

t  L i-i     l"^r_       2  J 


Here  n  is  the  generally  northerly  direction  parallel  to  the  satellite  sub- 
track .    A  typical  NIMBUS  III  subtrack  and  data  are  shown  in  Fig.   1.   . 

The  ballistic  D  values  and  values  of  u    and  v     calculated  from  all 

B  B 

available  rawinsonde  data  were  plotted  and  each  of  the  fields  analysed 
by  hand.    Values  were  then  interpolated  to  SIRS  subtrack.  positions. 
Since  the  points  with  SIRS  data  were  irregularly  spaced,  a  cubic  spline 
fitting  interpolation  routine  was  used  to  obtain  radiance  data  at  12  0  n.  mi 
intervals  along  the  satellite  subtrack.    Although  this  interpolation  pro- 
duced a  smoothing  of  the  radiance  values,  Wright  (1971)  has  shown  that 
the  standard  error  of  estimating  the  ballistic  D  value  was  changed  by  less 
than  two  meters  (about  four  percent) . 


<Q)   UNINTERPOLATED    SIRS-A 
DATA    POINTS 

©      INTERPOLATED   DATA     POINTS 


V_      COMPONENTS 
BN 


Fig.  1 ,  Typical  Nimbus  III  Satellite  Subtrack  Showing 
Uninterpolated  and  Interpolated  SIRS-A  Data, 
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B.     BALLISTIC  D  VALUE  APPROACH 

Correlation  of  the  ballistic  D  values  and  the  eight  cloud-corrected 

SIRS  radiances  led  to  the  regression  equation 

Dn  =  1.281    N.  -  5.243  N_  +  13.642  N     +  7.884  N     -  58.512  N 
B  1  2  3  4  5 

+  63.779  Nc  +  9.836  Nv  -  1.088  ND  -  1822.8  (13) 

6/8 

for  the  dependent  December  data.     For  the  January  data  the  corresponding 

equation  was 

D     =  -5.497  N,  -  58.457  N0  -  42.484  N0  +  32.944  N     -  91.296  Nr 
B  1  2  3  4  5 

+  40.878  Nr  -  10.980  N     +  14.980  N0  +  1792.3  (14) 

6  7  8 

Channels  2  through  6  contributed  significantly  to  the  D     in  the  January 

B 

sample,  while  the  leading  terms  for  the  December  sample  involved  only 
channels  5  and  6.     This  suggests  rather  poor  temporal  stability  and  the 
need  for  sequential  updating  of  the  dependent  sample.    As  discussed  by 
Elsberry  and  Martin  (1971),  the  earlier  period  was  marked  by  a  strato- 
spheric warming  with  a  generally  anomalous  tropospheric  flow  pattern  as 
well.    All  the  ballistic  winds  were  westerly  in  the  January  independent 
sample  (Table  1)  so  that  the  mean  V       was  equal  to  the  mean  of   |V       |. 
The  December  30  sample  contained  both  easterlies  and  westerlies  with 
an  average  magnitude  of  the  ballistic  wind  of  about  24  kts  .    As  a  result 
the  standard  deviation  of  the  December  sample  was  almost  30  kts  com- 
pared to  a  value  of  about  10  knots  in  the  more  uniform  January  sample. 

Estimates  of  Vr,>T  from  gradients  of  the  D  values  calculated  with  Eqs  . 
BN 

11 


w 
►J 

OQ 
H 


CD 
O 
C 
rd 


CO 

CD 
X3 
d 
O 

CO 

CJ 

£    fd 

fd   cc 

< 

s-t    Cl, 

M-l      (—1 

CD     Cn 
>     C 


cd 
Q 

CO 

C 


CO 

£> 

co 

G 

O 
-■-I 

>  g 

u  cj* 
.;;  w 

+-> 
co     C 

3   o 

CQ     w 

CD 

4-1     i-i 

O    CO 

o  ^ 
'd  o 

(0  £ 
a 

S  -a 
o  S 
u  * 


X 

o 

(0 

b 

s 

w 

S-i 

4-> 

o 

b 

w 

w 

to         _ 


CD     CD 

c   o 

jO     fd 
CL   "£ 

w  > 


E   cd 
o  o 

U  O 


c 

fd 
CD 


> 

CD 

Q 


CQ 


S3 
OQ 


CO    > 


OJ       CQ 

S    > 


ai 

^ 

CO 

CM 

CO 
■(-> 

CO 

+-> 

CO 

LO 

CO 

CO 

co 
O 

CO 


■ — i 


r-<        m 


cr> 

CO 

CN] 

CM 

CO 
■M 

CO 

4-> 

*4 

t^ 

CO 

CO 

i — i 

• 

• 

CT) 

vM 

CO 

+-> 

CO 
-M 

cn 

LO 

o 

LO  —l 


co 


CD 

co 


o 

CO 

CD 

CM 

CO 

CO 

cn 

LO 

oo 

CM 

CO 

cn 

LO 

CO 

CO 

PC 
W 
CQ 

2 

1 

W 

U 

OJ 

cn 

S3 

< 

o 

--H 

w 

• 

• 

i — . 

• 

• 

Q 

^r 

CO 

o 

r^ 

CN] 

CM 

"^ 

CO 

^ 


CO 


CM 
CO 

CO 

CO 

•4-> 

CO 

O 

-H 

O 

CO 

CO 
CM 


CM 


CD 


o 

fd 
o 
i-> 

a 
a 
< 


CD 

c 

3 

CD 

i—t 

-iH 

fO 

T3 

> 

rd 

u. 

Q 

o 

CD 

o 

fd 
> 


d 

CD 

fd 

a 


12 


(13)  and  (14)  using  SIRS-A  radiances  for  the  two  independent  samples  are 
also  shown  in  Table  1.    The  ballistic  wind  explained  variance  for  this  D 
value  approach  ranged  from  80  percent  for  the  December  sample  to  67  per- 
cent for  the  January  sample.    It  is  somewhat  surprising  that  the  standard 
error  of  the  latter  sample  was  only  5  kts  .  considering  the  magnitude  of 
the  ballistic  winds.     Such  errors  are  close  to  the  present  rawinsonde 
errors  in  high  wind  regions  .     Extrapolation  of  rawinsonde  winds  to  higher 
levels  for  integration  of  the  ballistic  wind  equation  can  also  introduce 
errors  of  this  size. 
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C.    GRADIENT  OF  N.  APPROACH 

1 

The  derivatives  X    =  — along  the  satellite  subtrack  were  corr- 

i     f     5  n 

elated  with  the  observed  values  of  the  normal  ballistic  wind  components 
V       .    For  the  December  dependent  data  the  resulting  regression  equation 
was 


VMT  =  -2.371    X,  -  11.456  X0  +  30.703  X0  -  1.548  X,  +  64.865  Xr 
BN  12  3  4  5 

-25.090XC    -  2.122  X.,    +2.251Xn    +24.42 
6  7  8 

Except  for  the  use  of  the  spline  fitting  technique  to  derive  radiances  at 
equally-spaced  latitude  points  in  the  gradient  calculation,  this  approach 
was  quite  direct.    As  a  result  of  differences  in  data-handling,  the  sample 
size  varied  slightly  for  this  method  compared  to  the  D  value  approach. 

Estimates  of  the  ballistic  winds  with  the  gradient  of  N.  approach 
(also  summarized  in  Table  1)  varied  considerably  for  the  two  independent 
samples.     Over  90  percent  of  the  variance  was  explained  with  a  standard 
error  of  about  eight  kts  for  the  December  sample.     Much  less  skill  was 
shown  for  the  results  in  January  since  only  about  40  percent  of  the  var- 
iance was  explained.     Even  though  the  standard  error  for  January  5  was 
only  about  10  kts,  this  is  relatively  large  compared  to  the  standard  devia- 
tion of  the  sample.    Evidently  the  larger  number  of  derivatives  required 
for  the  gradient  of  N.  approach  introduces  a  variability  in  the  results  with 
this  method.    Wright  (1971)  has  tested  this  method  without  interpolating 

the  radiances  .     In  both  the  gradient  of  N    and  the  ballistic  D  value 

i 

approaches,  the  use  of  variable  increments  caused  a  considerable 
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reduction  in  explained  variance  and  an  increase  in  the  standard  error.     In- 
herent smoothing  by  the  interpolation  scheme  seems  to  reduce  the  scatter 
in  the  radiance  data  introduced  by  the  independent  calculations  of  cloud 
correction  at  adjacent  subsatellite  scanspots. 
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IV.     SIRS-B  EXPERIMENT 
A.     DATA  AND  DATA  HANDLING 

The  SIRS-B  instrument  carried  onboard  NIMBUS  IV  was  noteworthy  for 
its  side-viewing  capability.     It  was  thus  anticipated  that  a  more  uniform 
distribution  of  radiance  data  would  be  available  to  permit  the  calculation 
of  both  horizontal  components  of  the  ballistic  wind.    A  30-day  period  in 
late  December  1970  and  January  1971  was  chosen,  as  the  high  winds  during 
this  period  would  provide  a  good  test  of  the  ballistic  wind  regression 
scheme.     Unfortunately  the  window  channel  was  inoperative  during  this 
period  so  only  the  seven  carbon  dioxide  channels  were  available. 

As  was  the  case  in  the  previous  experiments,  concurrent  rawinsonde 
data  and  corrected  clear  column  radiances  from  SIRS-B  were  required.    To 
permit  calculation  of  gradients,  the  data  were  chosen  from  areas  where 
several  SIRS  data  points  could  be  grouped  into  a  cluster.     These  clusters 
consisted  of  one  master  data  point  surrounded  by  three  to  six  other  "slave1 
points,  as  shown  in  Fig.  2.    In  areas  with  SIRS-B  coverage  the  normal 
spacing  between  data  points  was  450  to  650  km.    The  maximum  difference 
between  a  master  and  a  slave  was  arbitrarily  limited  to  1000  km.    As  only 
the  archived  NESS  clear-column  radiances  were  available,  the  data  cover- 
age was  generally  poor.    Data  were  restricted  to  the  latitude  belt  40N  to 
75N.    An  eight-day  period  was  used  to  collect  enough  data  to  derive  and 
test  regression  equations  analogous  to  those  in  the  SIRS-A  experiment.    A 

total  of  45  separate  data  clusters  comprising  194  points  were  found  in  the 
period  10-17  January  1971.     The  first  30  clusters  occurring  from  10-14 
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Fig.    2     Sample  SIRS-B  Data  Cluster 
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January  were  used  as  the  dependent  sample,  and  the  remaining  15  formed 
the  independent  sample. 

Calculation  of  both  the  east-west  (u)  and  north-south  (v)  components 
of  Eqs.   (11)  and  (12)  was  possible  with  the  SIRS-B  data  clusters.     In  Eq . 
(11)  the  gradient  of  the  ballistic  D  value  was  required  at  the  master  point, 
while  a  total  of  seven  gradients  of  the  radiances  in  the  cluster  was  necess- 
ary in  Eq .   (12).    Wright  (1971)  assumed  the  surface  between  each  master- 
slave  pair  was  linear  and  a  plane  was  fitted  through  the  data  of  each  pair. 
The  slope  of  each  plane  was  converted  to  an  east-west  and  north-south 
gradient,  and  the  average  of  such  gradients  for  all  master-slave  pairs 
(between  three  and  six  pairs)  was  taken  to  be  the  partial  derivative  of  the 
D  (or  N.)  value.    With  this  technique  the  required  partial  derivatives  could 
be  rapidly  calculated.    We  will  discuss  in  a  later  section  other  methods 
for  determining  the  D  value  gradient  at  the  center  of  an  irregularly-spaced 
data  cluster. 
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B.     BALLISTIC  D  VALUE  APPROACH 

Ballistic  D  values  were  calculated  using  radiosonde  temperatures 
from  stations  in  the  region  of  each  data  cluster  and  then  interpolated  to 
the  location  of  each  subsatellite  point  in  the  cluster.     Soundings  were 
extrapolated  only  to  1  mb,  rather  than  to  the  0.07  mb  level  used  in  the 
SIRS-A  study.     The  contribution  to  the  D  value  in  the  layer  between  these 
two  levels  was  generally  small,  and  always  based  on  extrapolated  data 
following  the  January  standard  atmosphere  lapse  rate. 

Correlation  of  rawinsonde-derived  D  values  with  the  radiance  data 

for  all  130  data  points  in  the  dependent  sample  produced  the  relation 

(15) 

D    =  -10.82  N0  +  15.23  N„     -47.41   N     +  95.96  N„  -30.45  N0  -1727.6 
B  ^  4  6  7  o 

This  equation  which  is  analogous  to  Eqs.   (13)  and  (14)  explained  only  81 

percent  of  the  variance  of  DD  in  the  dependent  sample  with  a  standard 

B 

error  of  120  meters.    The  large  standard  error  was  probably  due  to  the  long 
time  period,  as  well  as  to  large  variations  in  latitude  and  longitude.     It 
should  also  be  noted  that  the  window  channel  was  missing  during  this  per- 
iod and  only  five  of  the  remaining  seven  channels  gave  statistically  sig- 
nificant contributions . 

Radiances  from  the  15  clusters  in  the  independent  sample  were  used 
in  the  above  regression  equation  and  gradients  of  the  D  value  were  calcu- 
lated at  the  master  point  using  Wright's  technique.    These  SIRS-B  derived 
wind  components  were  correlated  with  the  rawinsonde  ballistic  winds 
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drawn  from  the  analysis  at  the  same  point.     The  mean  west  wind  compon- 
ent (see  Table  2)  of  21  kts  has  a  standard  deviation  of  19  kts  .    Again  the 
mean  of  the  absolute  value  of  this  component  shows  most  of  the  winds 
were  westerly.    Although  the  mean  north- south  component  was  near  zero, 
the  standard  deviation  of  v    was  about  15  kts.     In  the  D  value  approach 

D 

the  correlation  coefficient  for  the  north-south  component  was  just  slightly 
less  than  that  obtained  in  the  SIRS-A  results  for  the  normal  component. 
The  standard  error  of  the  estimate  of  9.5  kts  also  compares  favorably. 
However  the  estimate  using  the  D  value  approach  explains  only  27.6  per- 
cent of  the  variance  of  the  zonal  component.    A  standard  error  of  17  knots 
also  appears  excessive.     Some  of  this  error  is  due  to  the  oversimplifica- 
tion of  Wright's  technique  for  evaluating  the  gradient  from  the  irregularly 
spaced  array.    In  addition  the  distances  involved  in  the  gradient  calcu- 
lation were  large  compared  to  the  rather  closely-spaced  SIRS-A  data.    The 
larger  standard  error  of  the  D  value  estimate  from  Eq .   (15)  also  plays  a 
role. 
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C.    GRADIENT  OF  N.  APPROACH 

1 


As  shown  in  Eq .   (12)  the  u   and  v  components  of  the  ballistic  wind  may 
be  correlated  with  the  N.  gradients  evaluated  at  the  master  point  of  each 
cluster,  i.e.,  with 


X.  .        1        dNi  Y.  1      51 

i  =  -  —        - —         and         i    =   -r     t — 
f         dy  fox 


The  inclusion  of  the  Coriolis  parameter  (f)  reflects  the  normalization  of 
latitude  effects  in  the  sample.     Using  Wright's  technique  for  evaluating 
the  gradients,  equations  from  the  dependent  sample  (n  =  30)  were 


u     =152.2  X      -255.8  X     +  344.7  X     +  197.6  X      -695.7  X 

+  292.7  X     +  14.74 
8 

v     =  -133.4  Y      -121.0  Y     +  144.4  Y     -146.0  Y      +  568.8  Y 

-129.6  Y„    +  5595. 


Each  of  the  seven  channels   contributes  to  the  regression  in  this  method. 
For  the  dependent  samples  these  equations  explained  54  and  58  percent  of 
the  variance  in  the  u  and  v  components.    The  standard  error  for  each  com- 
ponent was  18  kts  . 

The  results  for  the  independent  sample  of  radiances  with  the  above 
equations  are  shown  in  Table  2.    Again  the  east-west  component  was  not 
estimated  as  well  as  the  north-south  component.    In  fact  the  results  for 

the  east-west  component  are  nearly  the  same  as  for  the  D  value  approach. 
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The  gradient  of  N.  approach  gave  slightly  poorer  estimates  of  the  north- 
south  component.    It  was  noted  in  the  SIRS-A  section  that  the  gradient  of 
N.  method  seems  to  be  consistently  less  useful  because  of  the  many  grad- 
ients to  be  calculated.    The  large  distances  over  which  the  gradients  were 
evaluated  probably  contributes  to  the  less  favorable  results. 

As  calculation  of  the  gradient  appeared  to  be  a  major  source  of  error 
with  the  SIRS-B  data  cluster,  in  the  next  section  we  examine  some  other 
methods  of  evaluating  the  gradient  at  the  master  point.    These  methods  are 
only  tested  with  the  ballistic  D  value  approach,  which  produced  superior 
results  in  both  the  SIRS-A  and  SIRS-B  experiments. 
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D.    ALTERNATIVE  METHODS  OF  DETERMINING  D  VALUE  GRADIENTS 

The  evaluation  of  the  ballistic  wind  components  at  the  master  location, 
corresponding  to  a  subsatellite  point,  is  dependent  on  the  east  and  west 
gradients  of  the  ballistic  D  values  through  the  geostrophic  relation.    Var- 
ious methods  can  be  used  to  calculate  these  gradients  given  the  D  values 
at  the  master  and  irregularly-spaced  slave  locations.    Wright  (1971)  used  a 
straight-forward  averaging  method  for  calculating  the  slopes  of  the  D  sur- 
face at  the  master  point.     More  general  relations  for  the  gradients  can  be 
obtained  by  expressing  the  set  of  D  values  in  the  region  of  master  point 

(X   ,  Y  )  as 
o       o 

2 
D  =  a.  4  a0  (x-x  )  +  a^  (y-y   )  +  a  .  (x-x  )  (y-y  )  +  a    (x-x  ) 
1         2  o  3  o  4  o  obo 

2 

+  ac(y-y  )    +  ...  (16) 

6         o 


Here  the  coefficient  a     represents  the  D  value  at  the  master  point  and  the 


and  the  gradients  of  D  may  be  written  as 

IT     =    a2  +  °4  (y-yo'  +  2  a5  (x-Xo' 

(17) 

W     =    a3+a4<X-Xo)+2a6(y-yo)  +  --- 

At  the  master  point  only  the  coefficients  a„  and  a     contribute,  but  their 

values  will  change  depending  on  the  number  of  terms  which  are  included 

in  the  expansion.    At  the  minimum,  there  must  be  one  master-slave  pair 
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for  each  additional  coefficient  added  in  the  expansion.    Thus  to  determine 
a     and  a     requires  at  least  three  master-slave  pairs.     If  more  than  this 
number  of  pairs  is  available,  the  coefficients  may  be  obtained  by  a  stat- 
istical regression  which  best  fits  (in  a  least-square  sense)  the  values  of 
D  at  all  the  available  points  in  a  cluster. 

The  multiple  regression  package  BIMED03R  of  Dixon  (1966)  was  used 
to  obtain  the  coefficients  a     and  a     with  three  degrees  of  approximations 
in  Eq  .   (16)    for  the  D  surface:    (1)  a  linear  regression  equation  including 

only  the  first  three  terms,    (2)  with  the  addition  of  one  second  order  term, 

2  2 

a.(x-x  )  (y-y  ),  and  (3)  with  the  addition  of  the  ar(x-x  )     and  a^(y-y  ) 
4         o  o  5         o  6         o 

terms,  but  not  the  a     term.    All  45  samples  used  by  Wright  (1971)  could  be 
used  with  the  first  (linear)  approximation,  but  only  33  samples  contained 
at  least  four  master-slave  pairs  and  could  be  used  with  the  second  approx- 
imation.    A  more  severe  reduction  in  sample  size  occurred  with  the  third 
approximation  as  only  nine  clusters  had  at  least  five  master-slave  pairs. 
This  reduction  in  sample  size  could  have  been  minimized  if  the  space 
(ten  degrees  latitude)  and  time  (less  than  two  hours  separation)  restric- 
tions set  by  Wright  (1971)  had  been  relaxed.     Nor  is  this  data  sample  typ- 
ical of  the  distribution  which  might  be  expected  if  the  technique  discussed 
here  were  to  be  used  operationally.    It  should  be  emphasized  that  the 
clear-column  SIRS-B  data  were  processed  for  a  different  purpose.    Many 
more  data  exist  on  the  SIRS-B  file,  but  the  selection  of  data  locations  and 

numbers  actually  processed  was  made  by  the  National  Environmental 

25 


Satellite  Service  for  calculating  temperature  profiles.    The  necessity  for 
a  different  processing  selection  for  use  with  this  technique  should  be 
obvious  . 

Since  at  the  master  point  the  coefficients  a     and  a     are  the  desired 
gradients  of  D  in  the  x-and  y-directions  ,  the  geostrophic  relation  may 
be  introduced  as 

b  D  fv  a_JD  fu 

d  x  g  3    y  g 

where  g  is  the  gravitational  acceleration  and  f  is  the  Coriolis  parameter 
at  the  master  location.     Both  g  and  f  are  constants  for  a  given  master- 
slave  cluster  and  thus  may  be  absorbed  in  the  regression  constant.    It 
might  also  be  noted  that  the  correction  for  longitudinal  increment  Ax  was 
important  since  the  data  sample  was  from  a  belt  between  40  N  and  75  N. 
Several  sources  of  error  may  be  distinguished  for  this  method.    The 
calculation  of  the  D  values  from  Eq.   (15)  using  the  satellite  radiances 
introduces  some  scatter  in  the  adjacent  values.    This  may  especially 
occur  if  very  different  corrections  to  the  observed  radiances  are  intro- 
duced due  to  different  cloud  cover.    Although  rather  large  standard  errors 
for  the  D  values  were  found  by  Wright  (1971),  often  the  majority  of  the 
points  in  a  cluster  would  have  the  same  sign  and  magnitude,  which  does 
not  affect  the  gradient.     Some  particular  examples  were  noted  near  the 
centers  of  deep  lows,  where  the  minimum  D  values  were  underestimated 
but  the  high  wind  speeds  were  correctly  calculated.    A  second  error 
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occurs  in  using  Eq  .   (16)  to  fit  a  surface  to  the  irregularly  spaced  data. 
In  a  few  cases  a  very  accurate  representation  of  the  D  values  calculated 
from  the  radiances  was  achieved,  but  the  corresponding  slopes  of  this 
surface  did  not  match  the  observed  wind  components  at  master  point. 
For  this  reason  an  optimum  surface  fitting  technique  should  also  be  able 
to  handle  D  values  from  Eq.   (15)  which  may  be  unrepresentative.     In  an 
operational  system  individual  clusters  should  not  be  treated,  but  rather 
an  objective  analysis  scheme  should  be  employed  so  that  the  calculation 
of  the  D  gradients  could  be  made  from  a  regularly  spaced  grid.     Unrepre- 
sentative D  values  would  be  filtered  by  the  objective  analysis  with  a 
good  first-guess  field  from  the  previous  analysis  and  proper  gross- 
error  checks  .    A  third  source  of  error  occurs  because  the  D  value  app- 
roach is  based  on  assumption  of  geostrophic  wind  components.    Thus  the 
wind  speed  will  be  overestimated  in  areas  of  cyclonic  curvature  since 
the  centripetal  acceleration  is  neglected. 

The  best  test  of  the  contribution  from  the  error  in  estimating  the  ball- 
istic D  value  from  the  satellite  radiances  should  be  to  calculate  the 
corresponding  geostrophic  wind  components  at  the  master  point  using  the 
observed  ballistic  D  values  at  all  the  points.     Thus  if  the  regression 

equation  (15)  for  D    was  perfect,  and  the  same  method  of  estimating  the 

B 

gradients  were  used,  the  ballistic  winds  calculated  using  the  satellite 
radiances  would  be  equal  to  those  using  the  radiosonde  ballistic  D  val- 
ues .     Because  the  radiances  are  from  a  single  instrument,  the  satellite 

27 


method  may  occasionally  produce  a  better  estimate  of  the  wind  in  regions 
of  sparse  or  inhomogeneous  rawinsonde  coverage.     In  the  following  tables 
the  ballistic  wind  estimates  from  both  the  satellite  radiance  data  and  the 
observed  radiosonde  data  will  be  compared  with  the  observed  ballistic 
wind  components  at  the  master  point.    As  noted  above,  three  degrees  of 
approximation  in  Eq.   (16)  will  be  used  in  testing  the  method  of  calculating 
the  ballistic  D  value  gradients  at  the  master  point,  which  is  related  to  the 
geostrophic  wind  component. 

The  first  test  (see  Table  3)  was  to  determine  the  ballistic  wind  compon- 
ents using  the  linear  regression  equation  approach  in  Eq .    (16).     Results  for 
both  the  dependent  (10-14  Jan  1971,  77=  30  samples)  and  the  independent 
(15-17  Jan  1971,  77=  15  samples)  cases  are  shown.    The  D  values  for  the 
latter  sample  are  derived  from  the  clear-column  radiances  using  Eq .   (15). 
Comparison  with  Wright's  slope-averaging  approach  (shown  in  parentheses 
in  Table  3)  shows  that  the  linear  regression  approach  raises  the  correlation 
coefficient  and  lowers  the  standard  error  of  the  estimate  in  the  cases  where 
these  were  least  satisfactory  using  the  Wright  approach.     However  the  re- 
verse is  generally  true  for  the  instances  in  which  the  Wright  approach 
showed  the  most  skill.    As  the  changes  are  small  there  appears  to  be  little 
statistical  difference  between  the  two  methods. 
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Table  3.    Correlation  coefficient  and  standard  error  of  the  estimate  for 
zonal  (u)  and  meridional  (v)  components  of  the  actual  ballistic  wind 
(kts)  calculated  by  a  linear  regression  equation  (see  text).     Corres- 
ponding values  obtained  by  a  slope-averaging  process  due  to  Wright' 
(1971)  are  shown  in  parentheses. 


Satellite 

u 
D  Values 


v 


Radiosonde 
D  Values 


Dependent 

Indepei 

ident 

V    = 

30 

V  = 

15 

.707 

16.8 

.558 

16.8 

(.661) 

(17.8) 

(.525) 

(17.0) 

.750 

16.6 

.779 

10.1 

(.775) 

(15.9) 

(.807) 

(9.5) 

.875 

11.5 

.796 

12.2 

(.860) 

(12.1) 

(.812) 

(11.8) 

.894 

11.3 

.576 

13.2 

(.875) 

(12.2) 

(.431) 

(14.6) 
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Comparison  of  the  satellite  D  value  approach  with  that  using  the 
radiosonde  D  values  shows  the  expected  decreased  accuracy  of  estimat- 
ing the  wind  speeds  for  all  samples  except  the  v  components  in  the  inde- 
pendent sample.    The  difference  in  accuracy  for  the  u  and  v  components 
appears  to  result  from  a  different  type  of  wind  regime  and  the  small  inde- 
pendent sample.    One  would  expect  there  would  be  no  preference  for 
wind  components  (as  shown  in  the  dependent  sample).    The  differences 
between  u  and  v  components  for  the  satellite-scan  and  rawinsonde- 
derived  values  certainly  implies  the  independent  sample  was  too  small. 
It  should  also  be  noted  that  use  of  the  independent  D  value  sample  from 
the  period  15-17  Jan  1971  was  a  rather  severe  test  of  the  temporal  stabil- 
ity based  on  a  dependent  sample  from  10-14  Jan  1971.    An  operational 
scheme  would  generally  be  updated  daily  using  the  most  recent  three-or 
four-day  period.     Finally  the  larger  degradation  of  the  zonal  component 
over  that  of  the  meriodional  component  suggests  the  importance  of  solv- 
ing for  D  regression  equations  in  latitudinal  bands. 

The  second  test  (Table  4)  in  determining  the  error  due  to  evaluating 

the  gradient  over  an  irregular  grid  was  to  include  one  second  order  term 

(x-x  )  (y-y  )  in  the  regression  equation.    Values  of  the  regression  coeff- 
o  o 

icient  and  standard  error  of  the  estimate  of  the  actual  ballistic  wind  are 
shown  in  Table  4  where  the  values  in  parentheses  are  for  a  homogeneous 
sample  from  the  linear  regression  technique.    For  this  small  sample 
rather  little  change  was  noted  with  the  inclusion  of  the  additional  term. 
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Only  three  of  the  eight  comparisons  showed  lower  standard  error.    It  is 
of  interest  to  note  that  the  apparent  improvement  over  Table  3  in  the  raw- 
insonde-derived  values  for  the  independent  sample  is  mainly  a  result  of 
choosing  a  different  sample. 
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Table  4  Correlation  coefficient  and  standard  error  of  the  estimate  for 
zonal  (u)  and  meridional  (v)  components  of  the  actual  ballistic  wind 
(kts)  calculated  with  a  regression  equation  including  an  xy  term. 
Corresponding  values  obtained  with  the  linear  regression  equation 
as  in  Table  3  are  shown  in  parentheses. 


Dependent  Independent 
N  =  21  N  =  12 

Satellite  .680        16.1  .549        18.3 

D  Values  U  1628)    (17.1)  (.588)     (17.7) 


v 


V 


.746        17.2  .757  9.9 

(.762)     (16.7)  (.772)      (9.7) 


Rawinsonde  .868       10.9  .888       10.1 

D  Values  U  (.903)     (9.4)  (.814)    (12.7) 


.882       12.1  .700       10.9 

(.869)     (12.8)  (.776)    (9.6) 
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For  the  third  test  the  (x-x  )(y-y  )  term  was  dropped  but  two  second 

o  o 

order  terms  were  substituted.     The  sample  size  (77  =  9)  in  this  case  is  so 

small  that  no  distinction  is  made  between  dependent  (4)  and  independent 

(5)  samples  in  Table  5.    As  noted  above,  the  inclusion  of  the    (x-x  ) 

o 

2 
(y-y  )  term  gives  mixed  results.     However  including  the  (x-x  )      and 
o  o 

.2 
(y-y  )     terms  does  lead  to  marked  improvement  in  each  case.     Since   the 
o 

sample  is  homogeneous  it  cannot  be  argued  that  these  results  merely  re- 
flect the  additional  data  available  in  these  nine  clusters.     More  samples 
would  be  required  to  establish  the  magnitude  of  this  improvement,  which 
is  consistent  with  the  better  representation  of  the  D  surface.    This  par- 
ticular sample  had  a  mean  zonal  wind  component  of  2  8.1  lets   with  a  stand- 
ard deviation  about  the  mean  of  22.7  kts  .    Corresponding  values  for  the 
meridional  wind  component  were  -2.7  kts  and  18.3  kts.    These  standard 

errors  of  the  estimate  of  less  than  10  kts  in  Table  5  for  the  regression 

2         2 
with  x     ,  y     term  thus  show  considerable  skill. 

One  auxiliary  product  of  the  BIMED03R  program  is  the  standard  error 

of  the  regression  coefficient,  which  in  this  case  represents  the  scatter 

in  the  u  or  v  components  .    That  is,  this  value  is  an  estimate  of  the  error 

made  in  not  being  able  to  fit  the  data  points  with  the  specified  surface 

(linear,  second  order,  etc.).     Large  values  thus  indicate  poor  confidence 

in  the  estimate  of  u  or  v  component.     In  the  large  majority  of  cases  for 

the  sample  shown  in  Table  5 ,  the  difference  between  the  actual  u  and  v 

and  the  estimated  u  and  v  components  was  smaller  than  the  standard 
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error  of  the  regression  coefficient.    However  it  must  again  be  cautioned 
that  the  D  values  may  be  fitted  very  accurately  by  the  specified  surface 
and  yet  not  predict  the  wind  components  if  some  of  the  D  values  are 
erroneous.    Thus  a  small  standard  error  of  the  regression  coefficient  is 
of  the  nature  of  a  necessary  but  not  a  sufficient  condition. 
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Table  5  Correlation  coefficient  and  standard  error  of  the  estimate  for 
zonal  (u)  and  meridional  (v)  components  of  the  actual  ballistic  wind 
(kts)  calculated  with  three  regression  equations  (see  text). 


2       2 
Equation  Linear  XY  Term  X    ,  Y    Term 


Satellite  u  .644  18.6  .647         18.5  .737  16.4 

D  Values  v  .785  12.1  .762  12.7  .894  8.8 

Rawinsonde  u  .890  11.1  .775  15.3  .960  6.8 

D  Values  v  .645  15.0  .735  13.3  .714  13.7 
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Ec     BALLISTIC  WINDS  USING  RAW  RADIANCES 

Examination  of  those  cases  with  poor  estimates  of  the  ballistic  wind, 
and  a  recent  paper  by  Hayden  (1971),  prompted  an  experiment  to  determine 
the  effect  of  the  cloud  correction  to  the  radiances.     In  several  cases  with- 
in the  independent  sample  the  cloud  corrected  radiances  in  the  channels 
sensing  the  upper  troposphere  and  stratosphere  varied  sporadically  about 
the  cluster  mean.    As  each  sounding  is  processed  independently,  the 
cloud  correction  could  represent  a  major  portion  of  the  gradient  between 
adjacent  soundings.     Hayden  (1971)  compared  SIRS-derived  soundings  with 
nearby  radionsondes  .    He  showed  the  quality  of  the  temperature  profiles 
from  the  SIRS-A  instrument  was  degraded  after  the  loss  of  the  714  cm 
CO„  channel,  which  suggested  a  similar  degradation  could  have  occurred 
with  the  loss  of  the  window  channel.    This  condition  may  have  been  one 
factor  in  the  less  accurate  ballistic  winds  obtained  with  the  SIRS-B  rad- 
iances.    Unfortunately  SIRS-B  data  with  the  window  channel  radiances 
were  not  available  to  directly  evaluate  the  effect  of  the  loss  of  the  win- 
dow channel. 

A  limited  experiment  with  the  D  value  approach  was  used  to  check 
the  effect  of  the  cloud  correction  in  the  SIRS-B  data  sample.    A  new  re- 
gression equation  was  calculated  relating  the  D  values  to  the  raw  rad- 
iances.   As  is  shown  in  Table  6,  the  explained  variance  in  the   depend- 
ent sample  was  0.769  versus  0.812  for  the  cloud-corrected  radiances. 
A  larger  standard  error  of  the  estimate  was  also  obtained  for  the  raw 
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Table  6.    Ballistic  parameters  calculated  using  raw  SIRS-B  radiances 
versus  clear-column  radiances  (shown  in  parenthesis). 


D  VALUE  CALCULATION 


Sample      Sample      Mean 
Size  D 


Depend-       145         -366.2       268.2 
ent 


Std         Correlation 

Explained 

Std 

Dev        Coefficient 

Variance 

Error 

268.2              .877 

.769 

131.7 

;269.8)           (.901) 

(.812) 

(119.2) 

Independ- 
ent 


78 


-330.4       165.8 


.765 


585 


107.5 


BALLISTIC  WIND  -  LINEAR  REGRESSION 


Depe 

mdent 
=  30 

Independent 
77=  15 

Component 

Correlation 
Coefficient 

Std 
Error 

Correlation 
Coefficient 

Std 
Error 

u 

.684 

17.4 

.600 

16.2 

(.707) 

(16.8) 

(.558) 

(16.8) 

V 

.535 

21.3 

.631 

12.53 

(.750) 

(16.6) 

(.779) 

(10. 1) 
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radiances.    The  test  with  the  independent  sample  resulted  in  a  reduced 
standard  error,  but  the  sample  standard  deviations  were  also  reduced. 

For  our  purposes  the  comparison  of  the  derived  wind  components  is 
the  more  important  criterion.    The  first  order  linear  regression  equation 
approach  discussed  in  the  previous  section  was  used  to  calculate  the 
D  value  gradient,  and  thus  the  derived  ballistic  wind  components.    Values 
in  parentheses  in  Table  6  were  abstracted  from  Table  3  for  comparison  with 
the  raw-radiance  winds.     Although  the  differences  in  the  correlation  co- 
efficient and  the  standard  error  were  rather  small  for  the  u  component  in 
the  dependent  sample,  those  for  the  v  component  were  much  worse.    The 
v  component  results  were  also  degraded  in  the  independent  sample,  while 
the  u  component  results  were  actually  improved.    It  is  clear  that  use  of 
the  cloud-corrected  radiances  is  important  for  the  v  component.     Not  much 
conclusive  can  be  stated  for  such  a  small  sample,  but  the  results  suggest 
that  some  of  the  reduced  accuracy  in  the  u  component  versus  the  v  com- 
ponent may  be  due  to  the  cloud  correction. 
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V.     CONCLUSIONS 

The  objective  of  this  research  was  to  demonstrate  the  feasibility  of 
calculating  ballistic  winds  using  SIRS  radiances,  and  to  explore   the  diff- 
erent approaches  for  calculating  the  ballistic  wind.    Ruggles  (see  section 
II)  has  suggested  the  mathematical  basis  for  a  direct  calculation  of  the 
ballistic  winds  by  evaluating  gradients  of  radiances  sensed  by  the  SIRS 
instrument.    In  a  second  approach  the  D  value  (deviation  from  the  stand- 
ard height)  times  the  ballistic  wind  weighting  factor  was  integrated 
through  a  deep  layer  of  the  atmosphere,  and  the  geostrophic  ballistic 
wind  was  calculated  from  the  gradients  of  these  ballistic  D  values. 

Only  components  normal  to  the  path  of  the  NIMBUS  III  satellite  could 
be  calculated  in  the  SIRS-A  experiment.    A  spline-fitting  routine  was  used 
to  interpolate  the  large  number  of  irregularly  spaced  "soundings"   to  equal 
increments  along  the  track.    Statistical  regression  between  the  ballistic 
wind  and  the  gradients  of  the  D  values  ,  or  of  the  radiances  ,  showed  very 
good  temporal  stability  when  tested  against  the  independent  sample  of 
30  December  1969  (see  Table  1)„    Although  the  regression  equations  with 
the  second,  smaller  independent  sample  on  5  January  1970  were  not  as 
stable,  the  test  suggests  standard  errors  of  less  than  10  kts  could  Le 
achieved. 

Both  ballistic  wind  components  can  be  calculated  from  the  "side-look- 
ing" SIRS-B  instrument,  but  the  method  of  calculating  the  gradient  re- 
quired clusters  of  SIRS  scan  points.    The  clusters  of  irregular, 
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widely-spaced  subsatellite  points  required  more  data  handling  and  re- 
stricted the  sample  size.    A  larger  standard  error  of  estimate  (see 
Table  2)  resulted  with  the  SIRS-B  data,  particularly  when  the  east-west 
component  of  the  wind  was  calculated .     Fitting  higher  order  approximations 
to  the  D  values  led  to  improved  estimates  of  the  wind  components  at  the 
central  point.    With  proper  data  handling  and  adequate  spacing  of  the 
processed  soundings,  it  would  appear  that  the  accuracy  should  be  com- 
parable to  that  for  the  SIRS -A  instrument. 

The  use  of  SIRS  radiances  for  deriving  temperature  profiles  over  a 
deep  layer  of  the  atmosphere  has  been  demonstrated  by  Smith  et  al 
(1970)  and  Hayden  (1971).    In  an  earlier  report  Elsberry  and  Martin  (1971) 
showed  that  the  ballistic  density  for  a  deep  layer  of  the  atmosphere 
could  be  calculated  directly  from  the  SIRS  radiances  .    The  capability  for 
viewing  a  deep  layer  of  the  atmosphere  is  the  basis  of  the  success  in 
directly  calculating  the  ballistic  wind.     Because  the  wind  depends  on 
gradients  of  the  integrated  density,  the  errors  in  the  wind  components 
are  proportionately  higher.     Nevertheless  the  proposed  technique  for 
determining  ballistic  winds  shows  considerable  promise,  especially  for 
the  oceanic  or  other  data-void  regions  . 
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clusters.   Because  the  geostrophic  ballistic  wind  is  proportional  to  a  hori- 
zontal gradient,  several  techniques  were  attempted  for  calculating  the 
gradient  at  the  central  point  of  the  irregularly  spaced  set  of  data.  With 
proper  data  handling  and  adequate  spacing  of  the  SIRS-B  soundings,  the 
accuracy  of  the  ballistic  wind  calculations  should  also  be  near  10  kts. 
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